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Abstract. In the early 17th century, Pierre de Fermat proposed the following problem: given three 
points in the plane, find a point such that the sum of its Euclidean distances to the three given points 
is minimal. This problem was solved by Evangelista Torricelli and was named the Fermat- Torricelli 
problem. A more general version of the Fermat- Torricelli problem asks for a point that minimizes 
the sum of the distances to a finite number of given points in R". This is one of the main problems 

"... 

in location science. In this paper, we revisit the Fermat-Torricelli problem from both theoretical 
and numerical viewpoints using some ingredients of of convex analysis and optimization. 



o 

(N 



O 



Key words. Fermat-Torricelli problem, Weiszfeld's algorithm, Kuhn's proof, convex analysis, sub- 
gradients 



AMS subject classifications. 49J52, 49J53, 90C31. 



u 
o 

S3 ■ 1 Introduction 



. The Fermat-Torricelli problem asks for a point that minimizes the sum of the distances 

!> . to three given points in the plane. This problem was proposed by Fermat and solved by 

fS| by the three given points is greater than or equal to 120°, the corresponding vertex is 

in 



the solution of the problem. Otherwise, the solution is a unique point inside of the triangle 
formed by the three points such that each side is seen at an angle of 120°. The first numerical 
CT} ■ algorithm for solving the general Fermat-Torricelli problem was introduced by Weiszfeld in 

1937 [18]. The assumptions that guarantee the convergence along with the proof were given 
by Kuhn in 1972. Kuhn also pointed out an example in which the Weiszfeld's algorithm fails 
^ , to converge; see [8]. The Fermat-Torricelli problem has attracted great attention from many 

researchers not only because of its mathematical beauty, but also because of its important 
applications to location science. Many generalized versions of the Fermat-Torricelli and 
several new algorithms have been introduced to deal with generalized Fermat-Torricelli 
problems as well as to improve the Weiszfeld's algorithm; see, e.g., [2, 9, 10, 12, 16, 17]. The 
problem has also been revisited several times from different viewpoints; see, e.g., [4, 5, 13, 19] 
and the references therein. In this picture, our goal is not to produce any new result, but 
to provide easy access to the problem from both theoretical and numerical aspects using 
some tools of convex analysis. These tools are presented in the paper by elementary proofs 
that are understandable for students with basic background in introduction to analysis. 

The paper is organized as follows. In Section 2, we proof the existence and uniqueness 
of the optimal solution. We also present the proofs of properties of the optimal solution 
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as well as its its construction using convex subdifferential. The advantage of using convex 
analysis when solving the Fermat-Torricelli problem has been observed in many books on 
convex and nonsmooth analysis; see, e.g., [3, 7, 15] and the references therein. Section 3 is 
devoted to revisiting Kuhn's proof of the convergence of the Weiszfeld's algorithm. In this 
section we follow the theme to proof the convergence given by Kuhn [8], but we include 
some ingredients from convex analysis to replace for some technical tools in order to make 
the proof more clear. 

Throughout the paper, IB denotes the closed unit ball of W 1 ; E3(x;r) denotes the closed 
ball with center x and radius r. 

2 Elements of Convex Analysis and Properties of Solutions 

In this section, we review important concepts of convex analysis to study the classical 
Fermat-Torricelli problem as well as the problem in the general form. We also present 
element proofs for some properties of optimal solutions of the problem. More details of 
convex analysis can be found, for instance, in [14]. 

Let || • || be the Euclidean norm in W 1 . Given a finite number of points ai for i = 1, . . . , m 
in W 1 , define 

m 

<p(x) = ^2 \\x - a,i\\. (2.1) 

8=1 

Then the mathematical model of the Fermat-Torricelli problem is 

minimize <p(x) subject to x € M n . (2-2) 

The weighted version of the problem can be formulated and treated by a similar way. 

Let / : M n — > M. be an real- valued function. The epigraph of / is a subset of W 1 x K defined 
by 

epi / := {(x,a) € M n+1 | x € M n and a > f(x)}. 
The function / is called convex if 

f(Xx + (1 - \)y) < Xf(x) + (1 - X)f(y) for all x,y € M" and AG (0,1). 

If this equality becomes strict for x ^ y, we say that / is strictly convex. We can prove that 
/ is a convex function on M n if and only if its epigraph is a convex set in M n+1 . 

It is clear that the function ip given by (2.1) is a convex function. 

Proposition 2.1 Let f : R n — > R be a convex function. Then f has a local minimum at x 
if and only if f has an absolute minimum at x. 

Proof: We only need to prove the implication since the converse is trivial. Suppose that / 
has a local minimum at x. Then there exists 5 > with 

f(u) > f(x) for all u G JB(x; 5). 
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For any x 6 K, one has that x^ = (1 — — )x + — x — > re. Thus, Xk € lB(x;5) when fe is 

k k 

sufficiently large. It follows that 

m < f( Xk ) < (i - + 

This implies 

and hence f(x) < f(x). Therefore, / has an absolute minimum at x. □ 

Proposition 2.2 The solution set of the Fermat-Torricelli problem (2.2) is nonempty. 

Proof: Let m := inf{<^(x) | x S R n }. Then m is a nonnegative real number. Let (x^.) be 
a sequence such that 

lim = m. 

k—^oo 

By definition, there exists ko £ N satisfying 

||xfc — ai | < (f(xk) < m + 1 for all k > &o- 

This implies < m+l+ Thus, (x&) is a bounded sequence, so it has a subsequence 
(xk e ) that converges to x € M n . Since ip is a continuous function, 

V?(x) = lim v?(xfc £ ) = m. 

i— >oo 

Therefore, x is an optimal solution of the problem. □ 
For two different points a, b € W 1 , the line containing a and 6 is the following set: 

C(a,b) = {ta + (l-t)b \ t € R}. 

Proposition 2.3 Suppose that ai for i = 1, . . . ,rn do not lie on the same line (not collinear). 
Then the function ip given by (2.1) is strictly convex and the Fermat-Torricelli problem (2.2) 
has a unique solution. 

Proof: Define (fi(x) := \\x — aj|| for i = 1, . . . ,m. Then <p = YllLi fi- For any x,y € R n 
and A € (0, 1), one has 

<Pi(\x + (1 - A)y) < A(^(x) + (1 - X)<Pi(y) for z = 1, . . . ,m. 

This implies 

^(Ax + (1 - \)y) < \<p(x) + (1 - \)<p(y). (2.3) 

On the contrary, suppose (p is not strictly convex. It means that there exist x, y £ R n with 
x ^ y and A £ (0, 1) for which (2.3) holds as equality. Then 

^i(Ax + (1 - X)y) = \ipi(x) + (1 - X)(fi(y) for i = 1, . . . , m. 
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Thus, 

\\X(x - at) + (1 - X)(y - ai )\\ = \\X(x - ai )\\ + ||(1 - X)(y - ai )\\ for i = 1, . . . ,m. 
Yix^di and y ^ aj, then there exists t% > such that 

tiA(x - a*) = (1 - A)(y - Oj). 

Thus, x — aj = 7i(y — Of), where ji := Since x ^ y, one has 7 7^ 1, and 

a,i = -t-^—x - —^—y G >C(x, y). 
1 — 7 1 — 7 

In the case where x = a, or y = a,}, it is obvious that € C(x,y). We have proved that 
aj € y) for i = 1, . . . , m, which is a contradiction. □ 

An element v £ M. n is called a subgradient of a convex function / : R n — >• M at x G M n if it 
satisfies the inequality 

/(a:) > f{x) + (v,x- x) for all x € M n , 

where (•, •) stands for the usual scalar product in M. n . The set of all subgradients of / at x 
is called the subdifferential of the function at x and is denoted by df(x). 

Directly from the definition, one has the following subdifferential Fermat rule: 

/ has an absolute minimum at x if and only if G df{x). (2-4) 

The proposition below shows that the subdifferential of a convex function at a given point 
reduces to the gradient at that point when the function is differentiable. 

Proposition 2.4 Suppose that f : M n — > M. is convex and Frechet differentiable at x. Then 

(Vf{x),x - x) < f{x) - f{x) for all x € R n . (2.5) 
Moreover, df{x) = {V/(x)}. 

Proof: Since / is Frechet differentiable at x, by definition, for any e > 0, there exists 5 > 
such that 

— e||x — x\\ < f(x) — f(x) — (Vf(x),x — x) < e\\x — x\\ whenever \\x — x\\ < 5. 

Define 

il>(x) = f(x) — f(x) — (Vf(x),x — x) + e\\x — x\\. 

Then ip{x) > ip{x) = for all x G S3(x;5). Since (p is a convex function, if)(x) > i^{x) for 
all x € M n . Thus, 

(Vf(x),x -x)< f(x) - f(x) + e\\x - x\\ for all x G R n . 
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Letting e — > 0, one obtains (2.5). 

Equality (2.5) implies that V/(x) G df(x). Take any t> G df(x), one has 

(v, x - s) < /(x) - /(x) for all x G R n 

The Frechet differentiability of / also implies that for any e > 0, there exists 5 > such 
that 

{v — V/(x), x — x) < e\\x — x|| whenever \\x — x\\ < 6. 

It follows that || v — V/(x)|| < e, which implies v = V/(x) since e > is arbitrary. Therefore, 
df(x) = {V/(x)}. □ 

The subdifferential formula for the norm function in the next example plays a crucial role 
in our subsequent analysis to solve the Fermat-Torricelli problem. 

Example 2.5 Let p(x) = \\x\\, the Euclidean norm function on M n . Then 

(B if x = 0, 

|^ | -jj — jj- 1 otherwise. 

x 

Since the function p is Frechet differentiable with Vp(x) = - — - for x 7^ 0, it suffices to 

ll x ll 

prove the formula for x = 0. By definition, an element v € W 1 is a subgradient of p at if 
and only if 

(v, x) = (v,x — 0) < p(x) — p(0) = \\x\\ for all x G ~R n . 

For i:=!)£ M n , one has < ||«||. This implies \\v\\ < 1 or v G iB. Moreover, if u G iB, 

by the Cauchy-Schwarz inequality, 

(v , x — 0) = (w, x) < ||u|| ||x|| < ||x|| = p(x) — p(0) for all x G M n . 

It follows that G dp(0). Therefore, dp(0) = IB. 

Solving the Fermat-Torricelli problem involves using the following simplified subdifferential 
rule for the sum of a nondifferentiable function and a differentiable function. A more general 
formula holds true when all of the functions involved are nondifferentiable. 

Proposition 2.6 Suppose that fi : M n — >• R for i = 1,2 are convex functions and fi is 
differentiable at x. Then 

d(/i + / 2 )(x) = <9A(x) + V/ 2 (x). (2.6) 
Proof: Fix any v G d{f\ + / 2 )(x). For any x G R n , one has 

(v,x- x) < fi(x) - fi(x) + f 2 (x) - / 2 (x) = /i(x) - fi(x) + (Vf 2 (x),x - x) + o(||x - x||). 
For any e > 0, there exists 5 > such that 

< (V/ 2 (x) — v, x — x) + /i(x) — /i(x) + e||x — x|| for x G iB(x; (5). 
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The convexity of f\ implies that this is true for all x. Letting e — > 0, one has 

< (V/ 2 (z) - v,x - x) + h(x) - h{x) for x € R n . 

By definition, v — V/ 2 (aO € dfi(x), and hence v € df\(x) + Vf 2 (x). We have proved the 
inclusion C. 

Since /1 = /1 + / 2 + (— / 2 ), by the previous proof, 

9/i (x) C + / 2 )(x) + V(-/ 2 )(S) = + f 2 )(x) - Vf 2 (x). 

This implies the opposite inclusion. □ 

Let us now use subgradients of the norm function to derive Torricelli's solution for the 
Fermat-Torricelli problem. Given two nonzero vectors u and v, define 

(u,v) 

cos{u,v) = - — ■ 
||li|| \\v\\ 

For x ^ a,i, let 

X-Oi . 

Vi = 7TZ 77, I = 1,2,3. 

\\x — CbiW 

Geometrically, Vi is the unit vector pointing in the direction from the vertex 04 to x. Observe 
that the classical Fermat-Torricelli problem always has a unique solution even if three given 
points are on the same line. In the latter case, the middle point is the solution of the 
problem. 



Proposition 2.7 Consider the Fermat-Torricelli problem given by three points 01,02,03. 

(i) Suppose x $l {01,02,03}. Then x is the solution of the problem if and only if 

cos(t>i,t> 2 ) = cos(u 2 ,^3) = cos(t> 3 ,-i>i) = -1/2. 

(ii) Suppose x € {ai, 02, 03}, say x = a\. Then x is the solution of the problem if and only 
if 

cos(v2,v-i) < -1/2. 

Proof: (i) In this case, the function ip given by (2.1) is Frechet differentiable at x. Since <p 
is convex, x is the solution of the Fermat-Torricelli problem if and only if 

V(p(x) = v 1 + v 2 + v 3 = 0. 

Since = 1 for i = 1,2,3, one has 

(vi,v 2 ) + (vi,v 3 ) 
(v 2 ,vi) + (v 2 ,v 3 ) 

{V3,vi) + (^3,^2) 

Solving this system of equations yields 

(vi,Vj) = cos(vi,Vj) = -1/2 for i / G {1,2,3}. 



= -1 
= -1 
= -1. 
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Moreover, if (vi,Vj) = — 1/2 for i ^ j, i,j £ {1,2,3}, then 

3 3 

IK + v 2 + -^3 II 2 = ^2 H^ll + K> u i) = o. 
i=i ij=i t ijtj 

It follows that v\ + v 2 + t>3 = 0. 

(ii) By the subdifferential Fermat rule (2.4) and the subdifferential sum rule (2.6), x = a\ 
is the solution of the Fermat-Torricelli problem if and only if 

G dipim) =]B + v 2 + v 3 . 

This is equivalent to |K + ^3 1 1 2 < 1 or |KI| 2 + IKI| 2 + 2(^2,^3) < 1. Since V2 and v 3 are 
unit vectors, we obtain 

(^2,^3) =cos(v 2 ,v 3 ) < -1/2. 
The proof is now complete. □ 

Example 2.8 Let us discuss the construction of the solution of the Fermat-Torricelli prob- 
lem in the plane. Consider the Fermat-Torricelli problem given by three points A, B, and 
C as in the figure above. If one of the angles of the triangle ABC is greater than or equal to 
120°, then the corresponding vertex is the solution of the problem by Proposition 2.7 (ii). 
Let us consider the case where non of the angles of the triangle is greater than or equal to 
120°. Construct two equilateral triangles ABD and ACE and let S be the intersection of 
DC and BE as in the figure. Two quadrilaterals ADBC and ABCE are convex, and hence 
S lies inside the triangle ABC. It is clear that two triangles DAC and BAE are congru- 
ent (SAS). A rotation of 60° about A maps the triangle DAC to the triangle BAE. The 
rotation maps CD to BE, so AD SB = 60°. Let T be the image of S though this rotation. 
Then T belongs to BE. It follows that /.AST = /ASE = 60°. Moreover, ZDS A = 60°, 
and hence ZBSA = 120°. It is now clear that ZASC = 120° and /BSC = 120°. By 
Proposition 2.7 (i), the point S is the solution of the problem. 



E 




3 The Weiszfeld's Algorithm 

In this section, we revisit Kuhn's proof [8] of the convergence of the Weiszfeld's algorithm 
[18] for solving the Fermat-Torricelli problem (2.2). With some additional ingredients of 
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convex analysis, we are able to provide a more clear picture of Kuhn's proof. Throughout 
this section, we assume that aj for i = 1, . . . , m are not collinear. 

The gradient of the function (p given by (2.1) is 

m 

Vip(x) = V] T l —, x £ {ai,a 2 , . . . ,a m }. 

^— ' \\x — a, 

i=l 

Solving the equation Vip(x) = gives 

Em a i 
i=1 \\ r — n-\\ 

:= b x). 



Em 
t=l 



\x — ai\\ 

For continuity, define F(x) := x for x G {ai, 02, ... , a m }. 

Weiszfeld introduce the following algorithm: choose a starting point Xo £ M n and define 

x fc+ i = F(x k ) for G IN. 

He also claimed that if x$ {oi, 02, • • ■ , a m }, where ai for i = 1, . . . , m are not collinear, 
then (xfc) converges to the unique optimal solution of the problem. A correct statement 
and the proof of the convergence were given by Kuhn in 1972. 

The proposition below guarantees that the function value decreases after each iteration; see 
[8, Subsection 3.1]. 

Proposition 3.1 If F(x) ^ x, then (p(F(x)) < <p(x). 

Proof: It is clear that x is not a vertex, since otherwise, F{x) = x. Moreover, F{x) is the 
unique minimizer of the following strictly convex function: 
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(*) = £ 



\z - ai|| 2 



U X - Oi\ 
1=1 



Since F(x) ^ x, one has g(F(x)) < g(x) = tp(x). 
Moreover, 

|2 

\x — ai " 



g(F(x))=J2 11 i_ n .; 
i=i 
in 

= E 



i=l 

x — OiW + \\F(x) — ad — x — a 



||x - aj| 



<p(x)+2((p(F(x))-(p(x)) + J2 



'"Fix) - aA\ - \\x - ai 



x - ai 
i=i 



It follows that 



W)) + £ (ll^-^ll-ll^-^ll) 2 < 2 ^(x). 

^— ' Hx — 1 



i=l 
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Therefore, tp(F(x)) < <p(x). □ 

The next two propositions show the behavior of the algorithm mapping F near a vertex 
and deal with the case where a vertex is the solution of the problem. Let us first present 
a necessary and sufficient condition for a vertex to be the optimal solution of the problem. 
It can be used to easily derive the result in [8, Subsection 2.1]. 

Define 

m 

Ui — a k 



u a i — a k\\ 
t=l,ij^k 

Proposition 3.2 The vertex a k is the optimal solution of the problem if and only if 

\\Rk\\ < 1. 

Proof: By the subdifferential Fermat rule (2.4) and the sub differential sum rule from 
Proposition 2.6, the vertex a k is the optimal solution of the problem if and only if 

Oedip(a k ) = -R k + B(Q,l). 

This is equivalent to < 1- □ 
Proposition 3.2 allows us to simplify the proof of the following result in [8, Subsection 3.2]. 

Proposition 3.3 Suppose that a k is not the optimal solution. Then there exists 5 > such 
that < \\x — a k \\ < 5 implies that there exists a positive integer s with 

\\F s (x) - a k \\ > 5 and \\F'- 1 (x) - a k \\ < 5. 
Proof: For any x, which is not a vertex, one has 

Em a i 
i=l 



F(x) 



\x — ai\ 



Em 
i=l 



Then 



F{x) - a k 



x — a,- 



\^m ®i ®k 

^ i= W* || x _ a-|| 



Em 
i=l 



X — Oj 



Thus, 



,. F(x)-a k ^i=i,i#* Hx-oill 

lim -r — = 1, [t- = K k . 

x^a-k \\x - a k \\ „ m ||x-a fc || 

1 + l^i=i,i^k 



\x — a; 



By Proposition 3.2, 

z-)-a fc II x — a k 



\\F(x) -a k \\ 
Inn — — = \\Ru\\ > 1. 
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Thus, there exist e > and S > such that 
\\F(x) — toll 

" ,, ' > (1 + e) whenever < \\x - a k \\ < 5. 

\\x - a k \\ 

and at ^ JB(a k , 5) for i ^ k. The conclusion then follows easily. □ 

We finally present Kuhn's statement and proof for the convergence of the Weiszfeld's algo- 
rithm; see [8, Subsection 3.4]. 

Theorem 3.4 Let (x k ) be the sequence formed by the Weiszfeld's algorithm. Suppose that 
x k {oi; 02, . . . , a m } for k > 0. Then (x k ) converges to the optimal solution x of the 
problem. 

Proof: In the case where x k = x k+ \ for some k = k$, one has that x k is a constant sequence 
for k > ko. Thus, it converges to Xk - Since F(xk ) = x ko and x ko is not a vertex, x ko is 
the solution of the problem. So we can assume that x k+ \ ^ x k for every k. By Proposition 
3.1, the sequence (ip(x k )) is nonnegative and decreasing, so it converges. It follows that 

lim (<p(x k ) - (p(x k+ i)) = 0. (3.7) 

k— >oo 

By definition, for k > 1, x k G co {a\, . . . ,a m }, which is a compact set. Then (x k ) has a 
convergent subsequence (xk e ) to a point z. It suffices to prove that z = x. By (3.7), 

lim (tp(x ke ) - ip(F(x kl )) = 0. 

I— >oo 

By the continuity, (f(z) = p(F(z)), which implies F(z) = z. If z is not a vertex, one has z 
is the solution of the problem, so z = x. Let us consider the case where z is a vertex, say 
ai. Suppose by contradiction that z/i. Choose 5 sufficiently small such that the property 
in Proposition 3.3 holds and IB{ai;5) does not contain x and an for i = 2, . . . , m. Since 
x k( — > a\ = z, we can assume without loss of generality that the sequence is contained in 

For x = x kl , choose qi such that x qi € IB{a\;5) and F(x qi ) £ ]B(a±;5). Choose an 
index kp > qi and apply Proposition 3.3, we find q<i > q\ such that x q2 € JB(ai;5) and 
F(x q2 ) £ JB{a\;5). Repeating this procedure, we find (x qe ) with x qe G JB(a\\5) and F(x qe ) 
is not in this ball. Extracting a further subsequence, we can assume that x qe — > q. By the 
procedure that has been used, one has F(q) = q. If q is not a vertex, then it is the solution, 
which is a contradiction because the solution x is not in JB(a\ \ 5). Thus, q is a vertex, which 
must be a\. Then 

r \\ F ( x q t ) ~ a lll 

hm — - — — = oo. 

£->oo \\X qe — ai || 

This is a contradiction according to Proposition 3.3. □ 

Acknowledgement. The author would like to thank Prof. Joel Shapiro for giving com- 
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